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1. Introduction 

The study of random planar graphs is a subject of interest in several areas of physics, 
ranging from two-dimensional quantum gravity, where planar graphs provide a natural 
discretization of space-time, to soft condensed matter where they can be used for instance 
to model fluid membranes. By random planar graphs, we mean here graphs embedded in 
the plane (or in the two-dimensional sphere), also referred to as fatgraphs in the physics 
literature or maps in mathematics. 

The statistical properties of random planar graphs have been studied by various tech- 
niques over the years. Beside the original combinatorial method due to Tutte [Q, and the 
quite technical matrix integral formulation P-Q, the most recent approach consists in a 
bijective enumeration 0, based on a one-to-one correspondence between planar graphs 
and suitably decorated trees, whose statistics is easily characterized by recursive relations. 
Apart from its simplicity, one of the advantages of the bijective approach is that it natu- 
rally keeps track of the distances between vertices on the graphs. This property was used 
in particular in Ref. to derive the statistical distribution for the so-called radius of quad- 
rangulations, and in Ref. to obtain the distribution of distances between two vertices in 
graphs made of even-sided polygons. It also gives access to a more refined study of the 
local environment of a vertex in random graphs, through the statistics of its number of 
nearest, next- nearest, ... neighbors. 

The purpose of this paper is precisely to use the bijective approach to derive the 
probability for a given vertex to have a number, say n^, of neighboring vertices and/or a 
number, say m^, of neighboring edges at a finite distance i. The case of nearest neighbors 
was already addressed in Ref. Q for the particular case of planar triangulations with no 
multiple edges and no loops. In this particular vertex with k neighbors may be 

replaced by a face with k sides, which makes the problem accessible to standard matrix 
integral techniques. This approach however breaks down for graphs with multiple edges 
or loops and is limited to the immediate neighbors. 

In this paper, we consider the case of both quadr angulations and Eulerian (i.e. face 
bi-colored) triangulations. In the former case, our analysis relies on a bijection described in 
Ref.|0] between rooted planar quadrangulations and so-called well-labeled trees, i.e. rooted 
planar trees with vertices labeled by non-negative integers, with a root vertex labeled 1 
and with the constraint that labels on adjacent vertices on the tree differ by at most 1. In 
the latter case, a similar bijection will be derived below between rooted planar Eulerian 
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triangulations and a more restricted class of well-labeled trees where the labels on adjacent 
vertices are required to differ by exactly 1. 

The paper is organized as follows. In Section 2, we present the correspondences 
between quadr angulations or Eulerian triangulations on one hand, and well-labeled trees 
on the other hand. In Section 3, we first recall known enumeration results for well-labeled 
trees (Sect. 3.1), and then show how to extract the statistics of vertices and/or edges at 
distance i from a given vertex by computing generating functions for suitably weighted well- 
labeled trees (Sect. 3. 2). When restricting to the case of the local environment (i.e. finite 
distances), we extend our study to the limiting case of infinitely large graphs (Sect. 3. 3). 
The following sections (4,5,6) are devoted to a detailed study of nearest neighbor statistics. 
In Section 4, we obtain algebraic equations for the associated generating function, which 
we use in Section 5 to calculate various probabilities and moments. Section 6 presents a 
heuristic interpretation of our results in the form of a phase diagram for a statistical model 
of embedded trees. We complete our study by presenting in Section 7 a shortcut to treat 
the general case of more distant neighbors. We gather a few concluding remarks in Section 
8. 

2. Geodesic distance in planar graphs: Prom graphs to well-labeled trees 

In this section, we first recall the construction of Ref. establishing a bijection be- 
tween on one hand rooted planar quadrangulations and on the other hand well-labeled 
trees, as defined above. This construction uses crucially the notion of geodesic distance 
between vertices in the quadrangulation and allows for an easy characterization of the 
number of neighbors of a given vertex. We then show how to adapt this construction 
to establish a bijection now between rooted planar Eulerian triangulations and a slightly 
different class of well-labeled trees. 

2.1. Quadrangulations 

Let us start with a rooted planar quadrangulation, i.e a planar map with all its faces 
of valence four (i.e. having four sides) and with a marked edge bordering the infinite face 
and oriented clockwise around the quadrangulation. The origin of the quadrangulation is 
the vertex from which the marked edge originates. For each vertex, we define its geodesic 
distance from the origin as the minimal number of edges of the paths linking the origin 
to this vertex. We may then label each vertex by its geodesic distance from the origin. 
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Fig.l: A sample rooted planar quadrangulation. An unambiguous represen- 
tation in the plane is obtained by having the marked edge adjacent to the 
external face and oriented clockwise around the graph. Each vertex is labeled 
by its geodesic distance from the origin of the marked edge (labeled 0). 

With this labeling, the origin is the only vertex labeled by 0, while all other vertices have 
positive labels. In particular, the nearest neighbors of the origin are the vertices labeled 
by 1 (see Fig.l). 

Each face of the quadrangulation may be classified, upon reading the four labels of 
the vertices around it in clockwise direction, in one of the two following categories (see 
Fig.2) : 

(i) "simple" faces with a sequence of labels of the form n-^n + l^n + 2^n + l 

(ii) "confluent" faces with a sequence of labels of the form + + 1 
Indeed, neighboring vertices have labels which differ by at most one by construction 

and the bipartite nature of the quadrangulation implies that the parity of the labels changes 
from any vertex to its neighbor. We now build a well-labeled tree with rules resulting from 
the above classification (see Fig.2): 

(i) For each simple face, we retain the edge n + 1 ^ n + 2 (in clockwise direction 0) 

(ii) For each confluent face, we draw a new edge linking the two vertices with labels n + 1 
The object resulting from the collection of these edges is a planar graph whose vertices 

are vertices of the original quadrangulation. Let us now show that it is indeed a well-labeled 

^ Here and throughout the paper, when we work in the planar representation, the rule must be 
reversed for the external face, i.e. the labels are to be read in counterclockwise direction around 
the graph. 
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Fig. 2: The two possible configurations of labels around a face in a quadr an- 
gulation: (i) a simple face and (ii) a confluent face. To each face is associated 
an edge (thick sohd line) as indicated below. The collection of these edges 
will eventually produce the well-labeled tree associated to the original quad- 
rangulation. 




(a) 



(b) 



Fig. 3: The impossibility of loops in the collections of edges as built from the 
rules of Fig. 2. Indeed, such a loop would have a vertex with minimal label n 
whose environment would be either (a) or (b). In both vertex of label 

n — 1 is present on both sides of the loop, which contradicts the definition of 
labels as geodesic distances from the origin. 



tree. Let denote the number of faces of the original quadrangulation, with therefore 
edges and N -\- 2 vertices. We note that N is also the number of edges of the graph 
resulting from the above construction as each face produces an edge and the resulting edges 
are all distinct (this is clear for confiuent faces and results from the orientation prescription 
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for simple faces). It is now sufficient to prove that the graph contains no foop. Suppose 
that such a foop exists, we pick a vertex with a minimal label, say n, along this loop. 
By examining all possible configurations (see Fig. 3), we see that this vertex is adjacent 
to at least one vertex labeled n — 1 on each side of the loop. This is a contradiction as 
the labels represent the geodesic distance: among the geodesic paths linking the origin to 
these vertices, one must cross the loop at a vertex with label less than n — 1, in conflict 
with the property that all labels on the loop are larger than n. The graph is therefore 
made of C > 1 connected planar trees, hence has N + C vertices. On the other hand, 
the number of vertices is at most + 1 as the origin cannot belong to a retained edge. 
We deduce that C — 1, hence the graph is a connected tree containing all the vertices of 
the quadrangulation but the origin. This tree is well-labeled as all labels are positive and 
labels on adjacent vertices diff'er by at most one by construction. 



Fig.4: The construction of the well-labeled tree associated with the rooted 
quadrangulation of Fig.l: the edges selected (a) according to the rules of Fig. 2 
produce the well-labeled tree (b). The tree is rooted at the vertex (labeled 
1) where the original marked edge ends. We have shaded the confluent faces 
in (a). 

The tree is naturally rooted according to the position of the marked edge in the original 
rooted quadrangulation (see Fig.4) 

So far we have associated to each rooted planar quadrangulation a rooted planar well- 
labeled tree. Let us now describe the inverse procedure. Starting from a well-labeled tree, 
we restore an origin by adding an extra vertex labeled 0. Non-intersecting edges are added 




(a) 



(b) 
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(a) 



(b) 




(c) (d) 

Fig.5: The inverse construction of Fig. 4 from the well-labeled tree (a) to the 
quadrangulation of Fig.l. We first restore the origin and connect it to each 
corner at all vertices labeled 1 (b). Within each of the newly created faces, 
further connections are made according to the rules explained in the text (c). 
The final step is to remove all edges connecting vertices with equal labels (d). 

between this vertex and each vertex labeled 1 with one connection per comer at this vertex 
(see Fig.5). The position of the root of the tree naturally translates into a marking of one 
of these newly created edges. The resulting map is made of a number of faces in which 
we apply the following construction independently. Reading the labels of the corners (as 
inherited from the label of the adjacent vertex) in clockwise direction around the face, we 
connect by a new edge each corner with label n > 2 to the first encountered corner with 



6 



label n — 1 , unless this latter corner is the immediate successor of the former and in such a 
way that the added edges do not cross each other. Finally, a planar rooted quadrangulation 
is obtained by erasing all the edges linking vertices with identical labels. We rely on Ref. 
for a formal proof that the above constructions are inverse of one-another. 

This bijection between rooted planar quadrangulations and well-labeled trees keeps 
track of the local environment of the origin. For instance, the vertices labeled 1 in the 
tree are the nearest neighbors of the origin, those labeled 2 the next-nearest neighbors 
and so on. This allows in particular to count the number of nearest neighbors of the 
origin as the number of vertices labeled 1 in the tree. Note that this number does not 
account for the multiplicity of the neighbors, which are all counted once even if they are 
connected by several edges to the origin. The number of nearest neighbors counted with 
their multiplicity is nothing but that of edges connected to the origin. In the language of 
well-labeled trees, this is the total number of corners labeled 1, or equivalently the sum of 
the valences of all vertices labeled 1. 

2.2. Eulerian triangulations 




Fig.6: A sample rooted planar Eulerian (face-bicolored) triangulation. An 
unambiguous representation in the plane is obtained by having the marked 
edge adjacent to the external face and oriented clockwise around the graph. 
This orientation induces an orientation of all the other edges as shown. Each 
vertex is labeled by its geodesic distance from the origin of the marked edge 
(labeled 0). This distance is now defined with paths respecting the orienta- 
tion. 
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We now intend to adapt the above constructions to the case of planar Eulerian trian- 
gulations, i.e. planar maps with all faces of valence three (hereafter referred to as triangles) 
and with an even number of triangles around each vertex (see Fig. 6). These triangulations 
are equivalently characterized by the property that their faces are bicolorable or that their 
vertices are tricolorable M. 





Fig. 7: The configuration of labels around a clockwise oriented (black) face in 
an Eulerian triangulation. To each black face is associated an edge (thick solid 
line) as indicated on the right. The collection of these edges will eventually 
produce the well-labeled tree associated to the original Eulerian triangulation. 



We start from a rooted such triangulation, with a marked edge bordering the infinite 
face and oriented clockwise around the triangulation. This fixes a natural orientation for 
every edge by requiring that orientations alternate around each vertex. This in turn gives 
a natural bicoloration of the faces in say black and white according to the clockwise or 
counterclockwise orientation of the adjacent edges. Taking as before for origin the vertex 
from which the marked edge originates, we define the geodesic distance of any vertex from 
the origin as the minimal number of edges of the paths linking the origin to that vertex and 
respecting the orientations. This defines a labeling of the vertices, such that the sequence 
of labels around all clockwise oriented (black) faces is of the form n—*n + l—*n + 2 (see 
Fig. 7). This follows from the fact that neighboring vertices have labels which differ by at 
most two while the tricolor ability of the vertices fixes the residue modulo 3 of all labels. 
We now build a well-labeled tree by a similar rule as for quadrangulations by retaining for 
each clockwise oriented (black) face with labels say + + 2 the edge linking 

the vertices labeled n + 1 and n + 2 (see Fig. 7). 

The resulting object is a well-labeled tree with the stronger constraint that adjacent 
vertices have labels differing by exactly 1. To see that, we follow the same line of reasoning 
as above. Let denote the number of black faces of the triangulations with therefore 2N 
faces, 3A^ edges and + 2 vertices. The graph resulting from the above construction has 



8 




Fig. 8: The impossibility of loops in the collections of edges as built from the 
rules of Fig.7. Indeed, such a loop would have a vertex with minimal label 
n which would be adjacent to a vertex labeled n — 1 on both sides of the 
loop, which contradicts the definition of labels as geodesic distances from the 
origin. 




(a) (b) 

Fig.9: The construction of the well-labeled tree associated with the rooted 
Eulerian triangulation of Fig. 6: the edges selected (a) according to the rules 
of Fig.7 build the well-labeled tree (b) 

N edges and contains no loop (see Fig. 8). It is therefore made of C connected trees and 
C = 1 since the total number of vertices AT -|- C is less than AT -|- 1 as the origin does not 

belong to any selected edge. 

The tree is naturally rooted according to the position of the marked edge in the original 
rooted Eulerian triangulation (see Fig.9) 

Let us now describe the inverse construction, which consists in iteratively building the 
black clockwise-oriented triangles. Starting from a rooted well-labeled tree with adjacent 
labels differing by exactly one, we construct the associated Eulerian triangulation step 
by step as follows. The step consists in adding an origin and linking it to the right of 
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(a) (b) 




(c) (d) 

Fig. 10: The inverse construction of Fig.9. Starting from the well-labeled 

tree, we first restore the origin and link it to the right of each 1 — > 2 edge (a). 
The links are then inflated into — > 1 2 black triangles (b). We proceed 
within each of the newly created white faces by linking the unique vertex 
labeled 1 within the face to the right of the 2 — > 3 edges (c). These links are 
inflated into 1 — > 2 — > 3 black triangles. The procedure is repeated until an 
Eulerian triangulation is obtained, which is already the case here (d). 

each edge with labels 1 and 2, oriented from 1 to 2. Each such link is inflated so as to 
create a clockwise oriented (black) triangle ^ 1 ^ 2. Beside these triangles, the original 
external face is split into a number of white faces (see Fig. 10). It is easily seen that the 
sequence of labels in clockwise orientation around any of these white faces is of the form 
0— >2— >2— >1 where all labels in-between the first and the last 2-labels are greater 
or equal to 2 and all increments but the first (0 — > 2) are ±1. In particular, if the sequence 
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is reduced to ^ 2 ^ 1, a white triangle has been created. We then proceed to the step 
1 which takes place independently inside each white face not reduced to a triangle. For a 
given such face F, each 2^3 edge having its right side adjacent to F is linked from this 
side to the unique vertex labeled 1 around F. As before, we inflate this link, thus creating 
a black triangle 1 — * 2 — ^ 3. As F is not reduced to a triangle, it is easy to see that at least 
one black triangle is created. Beside these additional black faces, the face F is split into 
smaller white faces. Among these faces, that containing the origin is a triangle ^ 2 ^ 1 
and the others of the form 1— s>3^---— s>3^2 such that all labels in-between the first 
and the last 3-labels are greater or equal to 3, and all increments but the first (1^3) are 
±1. These faces, if not reduced to triangles, have to be split again by repeating the above 
procedure. At step i, all faces are either black or white triangles or larger white faces of 
the form i^i + 2^---^i + 2^i + l where all labels in-between the first and the 
last i + 2-labels are greater or equal to z -|- 2 and all increments but the first {i ^ i + 2) 
are ±1. The procedure ends when all labels are exhausted and leaves us with an Eulerian 
triangulation rooted at the — 1 edge linking the origin to the root of the original tree. 

Let us now finally show that the two constructions above, respectively from an Eu- 
lerian triangulation to a well-labeled tree and conversely are inverse of one-another. It is 
straightforward that, starting from a well-labeled tree and building the associated triangu- 
lation, the latter is re-decomposed into the same labeled tree as the construction of black 
triangles always consists in adding edges of the form n —>■ n + 1 and n + 2 ^ ti to tree-edges 
of the form n + 1 —* n + 2. To complete the proof of the bijection it is now sufficient to 
show that the number of objects is the same within the two classes, namely on one hand 
rooted planar Eulerian triangulations with black faces, and on the other hand rooted 
well-labeled trees with N edges. The former number was first obtained in Ref. ]![] with the 
result 3 2^~^(2A^)!/(A^!(A^ + 2)!). The latter number will be derived below and is identical, 
which completes the proof. 

As for the case of quadrangulations, the bijection between Eulerian triangulations 
and well-labeled trees allows to keep track of the numbers of nearest, next nearest, etc ... 
neighbors of the origin in Eulerian triangulations. Note that the notion of neighbors now 
uses the edge orientations. As before, we have access to the number of neighbors both 
with and without their multiplicity. 
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3. Planar graphs viewed from a vertex: statistics of neighboring edges and 
vertices 



In this section, we use the above bijections to study the general statistics of edges and 
vertices at fixed geodesic distances from a given vertex (origin) in arbitrary quadrangula- 
tions and Eulerian triangulations. This is done via the enumeration of properly weighted 
well-labeled trees. We first recall how to perform this enumeration by use of recursion 
relations in the absence of weights, and then proceed to the weighted case. Finally, in 
the case of neighbors within a finite range of distances, we show how to derive algebraic 
equations for the various generating functions at hand. 

3.1. Enumeration of well-labeled trees 

The above bijections reduce the problem of enumeration of quadrangulations or Eu- 
lerian triangulation to that, much simpler, of well-labeled trees. This is done via the 
computation of their generating functions with a weight g per edge. More precisely, we in- 
troduce the generating functions Wj, j = 1, 2, ■ ■ ■, for rooted trees with positive labels and 
with a root labeled j. In view of the above bijections, the function Wi is indeed identified 
with the generating function for rooted quadrangulations (resp. Eulerian triangulations) 
with a weight g per face (resp. per black face). We are therefore only interested in Wi, 
but in practice, the scheme we propose involves the simultaneous computation of all Wj^s. 
From the constraint that adjacent vertices have labels differing by at most or exactly one, 
we deduce the following master equations: 

= l-giW,^,lw,+W,.,) (Quadrangulations) 

Wi = (Eulerian Triangulations) 

' l-g{Wj+i + Wj-i) 

obtained by inspection of the descendents of the root vertex (see Ref.|]T^). These recursion 
relations must be supplemented by the boundary condition Wq = and restricted to values 
of j = 1, 2, 3, ■ ■ ■. The solution is then unique as one readily sees by expanding eqs.(^T|) 
order by order in g and noting that all W^s are series of g with only non-negative powers 
by definition. 

The explicit solutions of eqs. ( p.l| ) were worked out in Ref. 0] , with the results 



W. = Jfz!f2±^ (Q.) 

Uj+iUj+2 

w. = Jfz!^^ (E.T.) 

Uj+iUj+3 
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(3.2) 



with 



1 - Jl - Via , , 

^ (3.3) 



and where 



«, = l-a:^ with = ^ = (1-12^7^ (Q.) ^34) 



x = \^, ^ = (l-8<7)4 (E.T.: 



1+ 

In particular, eqs.( p.2|) yield the following formulas for W\: 

4(1 + 2^ 

' 4(1 + 6*2)2 



(3.5) 



with Q as in eq.( p.4|) . For Eulerian triangulations, we easily extract from the above formula 
the number of well-labeled trees with edges 

where the subscript \gN refers to the coefficient of in the corresponding series. This 
number coincides with the number of rooted planar Eulerian triangulations with 2N faces, 
and allows to complete the proof of the bijection of Sect. 2, as announced. 

The function W of eq. (|3.3|) is the limit of Wj for large j. By a global shift of labels, 
the function Wj may be understood as counting trees with a root labeled 1 and with 
all labels > 2 — j. The function W therefore simply counts rooted trees with a root 
labeled 1 and with possibly negative integer labels. Interestingly enough, the function W 
may also be interpreted in the language of graphs. In the case of quadrangulations, it 
is easily identified as the generating function with a marked vertex and a marked edge, 
while Wi is the generating function with only a marked oriented edge (here we make 
use of the spherical representation of planar graphs for simplicity). If Qtv denotes the 
total number of quadrangulations on the sphere with faces (counted with their inverse 
symmetry factor), we have W\gN = {N + 2){2N)Qn while W^ilgW = 2{2N)Qn, hence 
WilgN = 2/{N + 2)W\gN = 23^cn/{N + 2) where cn = {^^)/{N + l) is the A^-th Catalan 
number, a result already obtained in Ref. [|ll| via tree conjugacy, and consistent with the 
solution (13.51), as checked in Ref. [fTl. 
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Similarly, in the case of Eulerian triangulations, the function W is identified with the 
generating function with a marked (origin) vertex and a marked oriented edge of type 
j —>■ j + 1 with respect to this origin and with a global orientation induced by that of the 
marked edge, while Wi is the generating function with a marked oriented edge (here again, 
we make use of the spherical representation of planar graphs for simplicity). If En denotes 
the total number of oriented Eulerian triangulations on the sphere with 2N faces (counted 
with their inverse symmetry factor), we have VF|gJv = {N + 2){2/3){3N)Eisf, as for a fixed 
origin, exactly 2/3 of the edges will be of type j ^ j + 1, while VFil^iv = {3N)En, hence 
WilgN = 3/{2{N + 2))W\gN =32^-^CN/iN + 2), a result already obtained in Ref.p via 



matrix models and in Ref. [12] via tree conjugacy, and also consistent with the solution 



3.2. Statistics of edges and vertices at fixed distances from an origin 

So far the counting functions have not kept track of the numbers of vertices and edges 
at a given distance from the origin in the graphs. This can be repaired by considering 
the same generating functions Wj as before but with extra weights, namely pi per vertex 
labeled i in the tree, and per edge adjacent to a vertex labeled i (an edge linking a vertex 
labeled z to a vertex labeled k thus receives the weight CKiCtfc). From now on, we shall use 
the shorthand notation Wj = Wj{{ai] pi}) for the generating function incorporating all 
the weights and pi, while the Wj's of previous section will be written as VFj({l;l}). 
Taking into account these new weights, the master equations become 

W- = (Q ) 

' l-gaj{aj+iWj+i + ajWj+aj-iWj-i) ^^^^ 

= z ^ r (E.T.) 

1 - gaj{aj+iWj+i + aj-iWj-i) 

which, together with Wq = 0, determine again the VF's completely as powers series in g. 

Let us now show how to extract from these generating functions the statistics of 
neighbors in quadrangulations of Eulerian triangulations. The neighbor statistics is entirely 
characterized by the following average: 



^iv({l;l}) 

ig,veg) 
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where the statistical sums extend over all pairs (Q, v) made of a quadrangulation Q (resp. 
an Eulerian triangulations Q) with N faces (resp. N black faces) on the sphere and a vertex 
V of Q. Here nj{Q,v) is the number of vertices at geodesic distance j from the vertex v 
while mj{Q,v) denotes the number of edges (j — linking a vertex at distance j — 1 
from V to a vertex at distance j from v. For instance, ni{Q,v) (resp. mi{Q,v)) denotes 
the number of nearest neighboring vertices (resp. edges) of v in Q. The factor \Aut{Q,v)\ 
is the usual symmetry factor of the marked graph on the sphere. 
The statistical sum A jv is related to Wi via 



«i^^Ar({«i;Pi}) = Wi{{aj] pj})\gr 



(3.9) 



The derivation with respect to ai amounts to counting the graphs Q (already marked at 
v) with a multiplicity 1711(0, v) equal to the number of edges originating from v. This 
corresponds precisely to counting graphs with a marked oriented edge on the sphere, i.e. 
rooted planar maps. 
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Fig.ll: The one-to-one correspondence between [j — 1, j) edges (represented 
in thick lines) in quadrangulations (top) or Eulerian triangulations (bottom) 
and edges adjacent to vertices labeled j in the corresponding well-labeled 
trees (dashed lines). 
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The relation (3^) follows from the fact that the total number of vertices at geodesic 
distance j from the origin is equal to that of vertices labeled j in the associated well-labeled 
trees, while the total number of edges {j — 1, j) in the graph is equal to the total number 
of edges adjacent to a vertex labeled j in the associated well-labeled tree. To see the 
latter property, note that the face adjacent to the {j — 1, j) edge and on its right contains 
exactly one edge of the tree adjacent to the vertex labeled j in this edge. This is done by 
inspecting all possible cases for quadrangulations and Eulerian triangulations (see Fig. 11). 
Conversely, to each edge adjacent to a vertex labeled j in the well-labeled tree, say of the 
form (j, i), we may associate a unique edge (j — l,j) of the graph. We have indeed the 
following possibilities for quadrangulations (respectively associated to the three situations 
of the first line of Fig. 11): 

(i) i = j — 1: the edge is the {j — 1, j) edge 

(ii) i = j: the edge (i, j) is a diagonal of a square (j — 1, j, j — l,i), and we pick the edge 
{j — 1, j) on its left 

(iii) i = j + 1: the edge (j, i) has a square on its right of the form {j — l,j,j + l,j), in 
which we pick the edge {j — 1, j) adjacent to (j, i) 

For Eulerian triangulations we only have two possibilities (respectively associated to the 
two situations of the second line of Fig. 11): 

(i) i = j — 1: the edge (z, j) is the {j — 1, j) edge 

(ii) i = j + 1: the edge (j, z) is adjacent to a black triangle on its right, of the form 
{j — 1, J, J + 1), in which we pick the {j — 1, j) edge 

This completes the bijection between the edges (j — 1, j) of the graph and those adjacent 
to the same vertex labeled j in the tree. 

In conclusion, inverting eq.(|3.9|), the knowledge of Wi as determined by eqs.( p.7| ) gives 
access to the statistical average F^r via 



N 



JV 



/o 1fW^i({«,i,---;i,i,---})l 

3.3. General solution for finite distances 



(3.10) 



From now on, we will concentrate on the study of the local environment of a given 
vertex. By this, we mean the statistics of edges and vertices within a finite range, say k, 
of geodesic distances from an origin. To this end, we may limit ourselves to a finite set of 
parameters {cti, ■ ■ ■ , ctfc; Pi, ■ ■ ■ Pfc} while all other ct's and p's are taken equal to one. Here 
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we wish to present a very general scheme for solving this problem in the form of algebraic 
equations for the generating functions. 
Introducing the functions 

Z, = a,W, (3.11) 
the equations (|3.7|) split into two sets of equations: a finite set 



Z. = fo^ 

' l-^7a,(^,+i + Z,+Z,_i) ^^-^ 

^ ^ ^ 



(3.12) 



valid for 1 < j < /c, and an infinite set 

1 



= (Q. 

' l-<?(Z,+i + Z,+Z,_i) 

= i-,(z,+, + z,_,) (^-^-^ 



(3.13) 



for all J > /c + 1. Remarkably enough, this latter infinite set of equations may be replaced 
by a single equation. Indeed, these equations may be solved order by order in g from 
the initial data Z^. Moreover, it is clear from the interpretation of the Zj = Wj for 
j > k + 1 that Zj — > as in ( p.3| ) when j — > oo. This is equivalently ensured by a single 
relation between Zk+i and Zk thanks to the existence of a discrete integral of motion 
f{Zj,Zj+i) = f{Zj+i,Zj+2) j > with 



f{x,y)=xy{l-g{x + y))-x-y (Q.) 

f{x, y) = xy{l - gx){l - gy) + gxy - x - y (E.T.) 

The convergence is ensured by writing 



(3.14) 



f{Zk,Zk+i) = f{W,W) (3.15) 



This equation is now supplemented by the finite set (|3.12|) of k equations, and upon 
elimination we may write a single algebraic equation for say Wi = Zi/ai. 

As an illustration, let us work out the case k = where all ct's and p's are one. We 
get a single equation /(O, Wi) = —Wi = f{W, W), which yields 

Wi = W- gW^ (Q.) 

(3.16) 

Wi = l + gW^{l-3gW) (E.T.) 
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in agreement with the results ( |3.5| ). 

The case k = 1 of the statistics of nearest neighboring edges and vertices wiU be 
studied extensively in the next section. 

For fixed /c, we have also access to the statistics of finitely distant edges and vertices 
from a vertex in infinite graphs, defined as the limit — > oo of the statistics at fixed A^. 
In particular, we may define 

T{{ai, ■ ■ ■ , ttfc; pi, ■ ■ ■ , pk}) = lim Ttv ({ai, ■ ■ ■ , afe, 1, ■ ■ ■ ; Pi, ■ ■ ■ , Pfc, 1, ■ ■ •}) (3.17) 

Jy —>oo 

Eq.( ^.lCi| ) translates into 

,s /o°' ^^^i({«,a2,---,afc,l,---;Pi,P2,---,Pfc,l,---})|smg 
r({ai, ■ ■ ■ , ttfc; pi, ■ ■ ■ , Pfcj) = -J— — 

(3.18) 

where the subscript |sing refers to the singular part of the corresponding function of g as g 
approaches the critical value gc = 1/12 (Q.) or gc = 1/8 (E.T.), proportional to {gc ~ g)^^"^ 
both in the numerator and denominator of eq. (|3.18| ), as we shall see in next section. 

4. Nearest neighbor statistics: general solution 

In this section as well as in the next two ones, we address the question of nearest 
neighbor statistics, for which explicit and particularly simple formulas may be written. 
Here we concentrate first on the case of quadrangulations and then repeat our analysis for 
Eulerian triangulations. 

4.1- Quadrangulations 



Combining the equations (|3.12|) and ( |3.15| ) for /c = 1 in the case of quadrangulations. 



we obtain the following third degree algebraic equation, with a = ai and p = pi 



{Wi-p){Wi{a-l)+p-l)-gaWi{W^a{a-l)+apWi+W{W-2))+2g'^aW^Wi = (4.1) 

We must pick the unique solution of eq.(^4.1|) with a g series expansion, and such that 
Wi = p + 0{g). We get 

Wi = p + gapil + ap) + g'^ap{2 + a + ap + 3a^p + 2ct^p^) + 0(5f^) (4.2) 
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Fig.l2: The rooted planar quadrangulations with up to two faces, and the 
corresponding weights ct per nearest neighboring edge (represented in thick 
sohd hne) of the origin (labeled 0), and p per nearest neighboring vertex 
(represented in grey, and labeled 1). These weights match the expansion 
(4.2) of VTi. 

up to order 2 in (7, in agreement with the direct enumeration of the quadrangulations with 
up to two faces (see Fig. 12). 

We may now study in detail the statistics of nearest neighbors in infinite quadrangu- 
lations by letting g ^ Qc = 1/12. Close to this critical point we write 



W 



1 + e 



(4.3) 



and substitute the expanded solution Wi = a+be+ce'^ +de^ +0{e^) into the cubic equation 
( ^T|) . At order in e, we find that 



F{a\a;p) = 36p(p-l)+4a(9-a-18p+9ap)+3a^(12(l-a)+a»+3aV(a-l) = (4.4) 



which fixes a as the unique solution with value 4/3ata = p=l, extracted from the 
solution ( |3.5| ) at 6* = ^/e = 0. Alternatively, this solution also satisfies a = p + 0{a) as a 
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series of a, and this latter property also fixes it entirely among the three possible solutions 
of F{a) = 0. Moreover we find at order one in e that 6 = 0, and the most singular 
(oc {qc — g)^^"^) term d, obtained at order 3 in e, is related to a via 

Wi{a; p)Ung = d = ^""^ (4.5) 

with F as in eq.( [4.4|) . Finally, let us perform the change of variable a —>■ A = aa in the 
integrals of eq.( p.l8|) . Using the explicit expression of a in terms of A: 

_ 3 12A + A^ + 24p-12-(2-A)v/(18-A)(2-A) 
2 36p{p + A-l)+A{3Ap + 3A^-4) 

as well as da/dA = —daF/{{A/a)daF — ad^F), we get 

(■A , 8(j;/a) / I \ 

Y-Tia-o) - ^y{y/c^)daF{y/c.\c.,p)-c.d^F{y/c.\c.-p) ^ 1 / - A \ 

^ ^ 2 IV 2- A i ^ ■ ^ 

Jo {ydaF{y\l-l)-d^F{y\l;l) \ / 

This remarkably simple formula is easily inverted into 

A - 2^Xr + 3)_ 

^-'(r + i)(r + 2) ^^-^^ 



itself mysteriously reminiscent of the form of the solution ( ^.21 ). To determine F, we finally 
note that a = A/a being a root of the cubic equation ( |4.4D , we may write a cubic equation 
for F as well: 

6F(1 + F)(3 + F) - a(2F(l + 4F + F^) + 3p(l + Tf{2 + F)) = (4.9) 

which fixes F uniquely by the condition that F = 1 for a = p = 1, or alternatively by F = 
when a = 0. 

4-2. Eulerian triangulations 



Combining the equations (|3.12| ) and (|3.15| ) for /c = 1 in the case of Eulerian triangu- 



lations, and taking again ai = a, pi = p, we get the following algebraic equation of degree 
three for Wi = Zi/a: 

{p - Wi){Wi{a - 1) + p - 1) + gaWi{W{W - 2) + (a - l)Wi{Wi + 1) 

(4.10) 

+ p(l + 2Wi - Wia - p)) + g'^aWiW'^il - 2W) + g'^aW^Wi = 
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Again, we must pick the unique solution of eq.( ^4.10| ) with a g series expansion, such 
that Wi = p + 0{g), whose first terms read 

Wi = p + gap + g'^ap{l + a + ap)+g^ap{3 + 2a + a^ + 2ap + 3a'^p + a'^p^)+0{g'^) (4.11) 

up to order 3 in g, in agreement with the graphs of Fig. 13. 
Close to the critical point we write 

g = -(l-e^), W = ^- (4.12) 

and substitute the expanded solution Wi = a+be+ce"^ +de^ +0{e'^) into the cubic equation 
(|]Tg) . At order in e, we find that 

F{a\a; p) = 32p{p - 1) + a(32 - 5a - 64p + 4ap{9 - p)) 

(4.13) 

+ 4a2((a - l)(a - 8) - ap(a - 2)) + 4a^a(a - 1) = 

which fixes a as the unique solution with value 5/4ata = p=l, extracted from the 
solution (|3.5|) at 6* = i/e = 0. As before, this solution also satisfies a = p + 0(a), which 
fixes it entirely among the three possible solutions of F{a) = 0. Again, we find at order 
one in e that 6 = 0, and the most singular (oc {gc — g)^^"^) term d is obtained at order 3 in 
e. It is related to a via 

with F as in eq.( [4.13|) . Finally, let us perform the change of variable a ^ A = aa in the 
integrals of eq.( p.l8|) . Using the explicit expression of a in terms of A: 



_ 9A-A' + 16p - 2Ap - 8 - (2 - A) ^(8 - A){2 - A) 
" 32p(p - 1) + 36Ap - AA^p - AAp^ + AA^-hA ^ ' ^ 

as well as da/dA = —daF/ {{A/a)daF — adaF), we get 

rA , 8{y/a) . — 

T = T(a- O) - y {y/^)daF{y/cc\c,,p)-c,d^F{y/c,\c,;p) _ 8- A ^ 

^ r^/^M, 8y \2-A ^ ' ' 

Jo {ydaF{y\l-l}-d^F(y\l;l) 

easily inverted into 

A = 2 , ^ J , (4.17) 

(r + i)(r + 3) ^ ^ 

also reminiscent of the form of the solution ( 3.2|) . We may again use eq. (|4.13| ) to write a 
cubic equation for F: 

4F(F + 2)(F + 4) - a(F + 1) (F(F + 5) + 2p(F + 2)(F + 3)) = (4.18) 

and F is fixed uniquely by the condition that F = 1 for a = p = 1, or alternatively by 
F = when a = 0. 
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Fig.13: The rooted planar Eulerian triangulations with up to 6 faces, and the 
corresponding weights a per nearest neighboring edge (represented in thick 
sohd hne) from the origin (labeled 0), and p per nearest neighboring vertex 
(represented in grey, and labeled 1). These weights match the expansion 
(4.11) of W^i. 
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5. Nearest neighbor statistics: probabilities and moments 

5.1. Quadrangulations 

We may extract from the solution of eg. ( [4 .91 ) various probabilities and moments for 
the numbers of nearest neighboring edges or vertices of a given vertex in infinite quadran- 
gulations. Taking for instance a = 1, the solution reads 

2 . V- . / 3 \ " f2n 



n>l 

in which we read the probability 

-(")^(il^)"e:) 

that a vertex have n nearest neighboring vertices. Analogously, taking p = 1, we get 



^, 1 /3(2 + a) \ 1 1 2 13 3 

in which we read the probabilities 1/3, 1/6, 13/108, ■ • • for having 1, 2, 3, ■ ■ ■ nearest neigh- 
boring edges. 

The formulae ( ^.1|) and ( p.3|) also yield the various moments of the distribution of the 
numbers ni and mi of nearest neighboring vertices and edges, upon writing p = (resp. 
a — e'") and expanding around u — (resp. v = 0). For instance we get 

W=3, (n?) = f, (n?) = ^ 

100 1372 (^•^) 

(mi) = 4, {ml) = — , {ml) = 

The value (mi) = 4 does not come as a surprise, as for each quadrangulation Q with 
faces, we have X]«ee "^1(^5 ''^) = 4A (twice the total number of edges), while there are 
N + 2 vertices, hence (mi) at = AN/{N 2) ^ 4 for large N. 

More generally, we have access to the joint probabilities for the numbers of nearest 
neighboring edges and vertices and to the associated mixed moments by expanding F as a 
double power series of a and p around (probabilities) or 1 (moments). 

Finally, we also have access to the conditional probabilities of having n nearest neigh- 
boring vertices given that there are no multiple nearest neighbors. This is achieved by 
taking p = t/a, in which case F = {a'^^~'^^t'^^), and then sending ct to retain only 
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the configurations where mi = ni (as mi > ni always). This results in the generating 
function n(t)/n(l) for the above conditional probabilities, where 

, t , if 118 -t \ 1 7 n 53 o , , 

U{t) = lim r(a; -) = - \ 3 = -t + —t^ + —t^ + ■■■ (5.5) 



' a' 2 \ \ 2-t / 3 54 972 



and in particular we get the probability of having no multiple nearest neighbor 11(1) = 
(\/l7 — 3)/2. We note that from eq.( [4.7| ) we have 11(74) = F. This relation is not surprising, 
once understood in the language of well-labeled trees. Indeed, the trees contributing to 
n(/:) are those associated with graphs without multiple nearest neighbors of the origin, 
i.e. those in which the vertices labeled 1 are all terminal (univalent) vertices and receive 
a weight t. Replacing this weight thj Zi = aWi reconstructs the full generating function 
contributing to F. For large graphs, this corresponds to the substitution t A in II. 

5.2. Eulerian triangulations 

Repeating the various specializations of previous section of the solution to eq.( [4.18|) , 
we get for a = 1: 



1 f j27-2p ^\ 4 _ 56 2 848 3 



leading to the probabilities 4/9, 56/243, 848/6561, ■ ■ ■ that a vertex have 1, 2, 3, ■ ■ ■ nearest 
neighboring vertices. For p = 1, we get: 

F(a; 1) = ^=i= - 1 (5.7) 
V 4 — 6a 

in which we read the probability 



Pim) = [ (^'^] (5.8) 



16/ \ m 

for having m nearest neighboring edges. Remarkably, these probabilities for nearest neigh- 
boring edges in Eulerian triangulations coincide with those of eq.( p.2| ) for nearest neigh- 
boring vertices in quadrangulations. Note also that the probability ( p.8| ) is very similar to 
that obtained in Ref. for ordinary triangulations with no multiple edge and no loop. 
The corresponding moments read 

12 , 2v 1212 , 190236 
^ ' 5 ' ^ 125 ' ^ 3125 

(mi) = 3, {ml) = y, (mf) = — 
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The value (mi) = 3 can be obtained by explicitly computing {mi)^ = 3N/{N + 2) ^ 3 
for large A^. 

Finally, we also have the conditional probabilities of having n nearest neighboring 
vertices given that there are no multiple nearest neighbors. Repeating the limiting process 
of previous section we arrive at 

^ , t^ /8-t „ 3 39 2 267 o 

U(t) = m r a; - = W 2 = -t H H H 5.10 

^ ' c^o ^ ' a' V 2-t 8 256 4096 ^ ' 

and in particular we get the probability of having no multiple nearest neighbor 11(1) = 
77-2. As before, we have U{A) = T of eq.(P^). 

6. Heuristic interpretation: Phase diagram 

The results of previous section display a number of critical values of a and p at which 
the various generating functions become singular. For quadrangulations we have singular 
points at (a = l;p = 4/3) and {a = 6/5; p = 1), easily read off eqs.(pTT|) and ( |5.3| ). The 
value p~^ = 3/4 (resp. = 5/6) governs the exponential decay of the corresponding 
probabilities, behaving as (3/4)"^ (resp. as (5/6)"^). Similarly, for Eulerian triangulations, 
we have the particular points {a = 1; p = 3/2) and (a = 4/3; p = 1), read off eqs.(|5^) 
and ( p.7| ). These points are particular points on a more general critical curve in the (a, p) 
plane, whose location may be easily obtained from the following heuristic argument. 

Labeled trees may be thought of as trees embedded in a one-dimensional target space, 
by viewing the labels as the positions of the vertices in the target space. In this language, 
the generating function Wj corresponds to trees with a root at position j. Without the 
constraint that the labels remain positive, all Wj's are equal to W of eq. (|3.3|) , with W\gN = 
3^CAr (Q.) or 2^cn (E.T.), counting arbitrary rooted planar trees with edges (in number 
cn) together with 3^ (resp. 2^) embeddings accounting for the 3 (resp. 2) possibilities of 
labeling a vertex according to the label of its immediate ascendent vertex. 

The positivity of all labels simply restricts the target space to a half-line, by introduc- 
ing a hard wall at the origin, whose effect is to reduce the entropy when compared with the 
no-wall situation. On the other hand, the weights a and p of Sect. 4 act as energy weights 
when an edge or a vertex touches the wall, which may help restore the balance between 
the reduction of entropy and a possible gain in energy. The critical curve corresponds to 
an unbinding transition at which the balance is realized. Its location may be obtained 
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Fig. 14: Phase diagram in tlie (p, a) plane for rooted trees embedded in a one- 
dimensional target space with a wall and weights a (resp. p) per edge (resp. 
vertex) touching the wall. A critical curve separates a region where the tree 
is attracted by the wall from a region where it is repulsed by it. This curve 
possesses a higher order (tricritical) point represented by a black dot. The 
corresponding leading singular behavior of Wi is indicated. The indicated 
vahies of a and p refer to well-labeled trees associated to quadrangulations 
(Eulerian triangulations) . 

as follows. For each vertex of the tree reaching position 1, we have an energy gain pa^ , 
where k is the number of nearest neighbors of the vertex on the tree. On the other hand, 
the entropy loss is (2/3)^^"^ (Q.) or l/2^~^ (E-T.) as the k — 1 descendents of this vertex 
have lost one of their 3 (resp. 2) possible positions. In an infinite rooted planar tree, the 
probability of having k neighbors is 1/2 hence the balance is reached when 



This displays the critical lines 2(3 — a) — 3pQ; (Q.) or 4 — o; = 2pa (E.T.) separating a 
phase in which the tree is "attracted" by the wall from a phase where it is "repulsed" by it 
(see Fig. 14). All pairs (a; p) mentioned above belong to these critical hues. A special point 
on these lines corresponds to the exact balance for each separate value of the valence A;, 
namely when pa^ = (3/2)'=-i (Q.) or pa^ = 2^'^ (E.T.) for aU /c > 1, i.e. a = 3/2, p = 2/3 
or CK = 2, p = 1/2 (E.T.). At this special point, the effect of the wall is exactly suppressed 
and the half- line solution coincides with the no- wall solution, namely Zj = W for all j > I. 




(Q-) 



(6.1) 




(E.T.) 
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Therefore at this special point, the singularity of Wi = pW is proportional to {qc — qY^'^, 
instead of the generic singularity oc {qc — g)^^"^ already observed away from the critical 
line. It is easily checked that the generic singularity along the critical line is of the form 
{gc — g)^^^: easily obtained by expanding the equations ( ^71] ) and ( [4.1C1|) now in powers of 
9 = i/e and solving order by order for the coefficients of Wi = A + B6 + C6^ + D6^ + ■ ■ ■ 
These various critical behaviors translate into different asymptotics for the weighted sums 
over graphs of size A^, of the form g~^ /N'^ with 7 = 5/2 (most generic case, away from 
the critical line), 7 = 7/4 (generic on the critical line) and 7 = 3/2 (special point on the 
critical line). 

Beside the above heuristic argument, the critical lines may also be derived by imposing 
that both F{a\a;p) = and daF{a\a; p) = (with F as in eqs.( |4.4[ ) or ( [4.13| )), while the 
special (tricritical) point corresponds to having in addition d'^F{a\a; p) = 0. 



7. Generalization to more distant neighbors 

7.1. General scheme for getting F 

The results of Sects. 4 and 5 may be generalized so as to include more distant neigh- 
bors, say up to k-th neighboring edges and vertices. This involves explicitly solving the 



equations (|3.12| ) and (|3.15| ). The resulting algebraic equation for Wi however grows in 
degree with k, and one still has to integrate once with respect to ai to finally obtain 
the statistical average F of eq. ( |3.18| ). Nevertheless, we have found a particularly simple 
alternative way of getting F by relating it to a single quantity t parametrizing all Zj's at 
the critical point. More precisely, when g = g^ with gc = 1/12 (Q.) or gc = 1/8 (E.T.) 
and W = Wc = 2, x = Xc = I (as readily seen from eqs.(^.3D and (pi^)), the convergence 
condition ( |3.15D is trivialized by a change of variables on = borrowed from the exact 
solution. Namely, if we write 



_ o {t+j){t+j+3) fp^ 

Z^^Qit,), With (7.1) 

^(t+j + l)(t+j+3) \^-^-) 

for all j > k, then the condition f{Zj, -^?_|_i) = /(2, 2) amounts to tj = tj+i, hence tj = t 




for all j > k. The value of t is obtained via the set of k equations (|3.12| ), which we feed with 
the two "final conditions" Z^ = (kii) and Z^_^-^ = (k+i{t), allowing to express backward 
all Zj^s with j < k as rational functions of t and to finally get one algebraic equation for t 
with ck's and p's as parameters. The proper solution is then uniquely characterized by the 
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fact that when aU a's and p's are one, we have t = 0, as can be read off the solution 
at criticahty, where W = 2 and Uj oc j. 

Quite magicaUy, it turns out that t may be identified with the desired function F, up 
to a translation independent of ai. More precisely, we have 

r({ai; Pi}) = 02, ■ ■ ■ ; Pi, P2, ■ ■ ■}) - t({0, ^2, ■ ■ ■ ; Pi, P2, ■ ■ ■}) (7.2) 

as r must vanish when ai = (note that when cti = 0, Zi = and t does not 
depend on pi either). This remarkable identity may be checked at /c = 1 where 
t({0, 1, 1, ■ ■ ■ ; pi, 1, 1, ■ ■ ■}) = t({0, 1, 1, ■ ■ ■ ; 0, 1, 1, ■ ■ ■}) = —1 (the solution is the same as 
when all ct's and p's are one except for a shift in the indices) and the relation F = t + 1 is 
read off eqs. (|4.8|) and ( [4.17|) with A = . For arbitrary k, the identification ( [7.2|) will be 



proved in the next section thanks to scaling arguments. 

The value of F may therefore be extracted as follows: let G{t) = G{t\{ai; pi}) = 
denote the algebraic equation satisfied by t. Then if we take cti = we get an algebraic 
equation G'o(to) = for the translation parameter to = ^({0, 02, ■ ■ ■ ; Pi, P2, ■ • ■}) of eq.([7.2|). 
We therefore simply have to eliminate to between the two equations ^(F + to) =0 and 
Go{to) = to get an equation for F. This equation might be sometimes factorized into 
smaller degree equations, but the initial condition F = 0{ai) always selects only one of 
them. For illustration, we will recover in Sect. 7.3 below the results for F at /c = 1, and 
further present the solution for /c = 2. 

7.2. The reason why t = T + to 

In this section we show that the parameter t occurring in the changes of variables 

(t + J + l)(t + J + 2) ^^-^ 

^ {t + j + l){t + j + 3) ^ ' 

coincides with the statistical average F up to an additive constant independent of cti. We 
simply have to show that dt/dai = dV/dai. Let us first find a simpler characterization of 
t, by considering Zj at large j. Expanding Zj at large j > k, we find 



z] = ^'^ + ^it + l)+0{j,) (Q.) 

Z] = 2-^ + ^{t + 2) + 0{-J (E.T.) 
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(7.4) 



therefore t enters only the subleading correction oc of Zj. Let us now go back to 
the combinatorial interpretation of Zj (for large enough j) as a series of g, before taking 
g gc- We see that dZj/dai is the generating function for well-labeled trees with a root 
at position j and a marked oriented edge originating from a vertex at position 1. In terms 
of graphsU, consider the generating function Gj({a^; pi}) for graphs with a marked vertex 
(origin) and a marked edge from a vertex at geodesic distance j — 1 from the origin to a 
vertex at geodesic distance j from the origin. Then we have 

Gj{{ai, ^2, ■ ■ ■ , ttfc; Pi, P2, ■ ■ ■ , pk}) 

(7.5) 

= ^j({ai,a2, ■ ■ ■ ,afc;Pi,P2, ■ ■ ■,Pfc}) - ^j-i({a2,a3, ■ ■ ■ ,afc;P2,P3, ■ ■ • ,P/c}) 

Indeed, in the tree language, Gj generates the well-labeled trees with root at position j 
and which do reach the position 1, therefore expressed via ( [7.5| ) by subtracting from the 
generating function Zj{{ai;pi}) the contribution of trees which do not reach position 1, 
also equal to that of trees with root at position j — 1, provided we shift the whole picture 
by —1, thus removing the terms pertaining to position 1. Taking a derivative of 



eq.(|7.5|) with respect to ai, we see that dZj/dai = dCj/dai. With this identification, we 
have 



When ^ oo, we deduce that 



dZj 



aid^lsing _ {mimjllaT'p7!)_ ^7 



W^{{l■,l})Uns (mi) 

where (■ ■ ■) with no index refers to the large limit. For large j we expect the average 
in the numerator to factor into {mj){mi Y[ c^T^ pT) weights are localized around the 

origin {k is kept finite). We arrive at the relation 

hm ^'""^ = ^i({l;l})lBmg dr _ 2 dV ^^g^ 
j^oo {rrij) {mi) dai 3 dai 

for both quadr angulations and Eulerian triangulations: indeed we have used (mi) = 4 and 
M^i({l;l})|,ing = 8/3 (Q.) while (mi) = 3 and M^i({l; l})|sing = 2 (E.T.) from eqs.(P) 
and ( [4.141 ) at a = p = 1. 



^ Here and throughout the appendix, the mention graph refers to either quadrangulations or 
Eulerian triangulations. 
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Both {rrij) and dZ j / dai\s\ng may be computed by considering the usual scahng hmit 
reached by sending simultaneously g Qc and j — > oo, say by setting 

g = g,il-6^), j = ^ (7.9) 

and letting 9^0 (here plays the role of a correlation length for the models at hand). 
In this limit, we have 

Zj = 2(1 - e^U{r) - e^V{r) ) (7.10) 

Upon substituting this expression for Zj into eqs.( |3.1^ ), we obtain at order 4 and 5 in 6* 
the following differential equations 

U" = 3 (W^ - 1) V" = 6UV (Q.) 

(7.11) 

U" = 2{U^ -1) V" = AUV (E.T.) 

The behavior of the scaling functions U and V when r — corresponds to taking g = gc 
at fixed j, hence we read off the expansions (^]3) °f following leading behaviors 

U{r) ~ 2/r2 and V(r) = -4{t + D/r^ (Q.) while U{r) ~ 3/r2 and V(r) = -6{t + 2)/r^ 
(E.T.). The functions U and V are further fixed by noting that Zj W = 2/(1 + 6'^) = 
2(1 — 9"^ + 0{9'^)) at large j, hence U{r) — > 1 and V(r) — when r oo. This results in 





while 



W(r) = + 1 (Q.) 

sinh lA/lr 



W(r) = —-2— + 1 (E.T.) 
smh [r) 

V(r)= (t + 0W'(r) (Q.) 

V(r) = (t + 2)W'(r) (E.T.) 



(7.12) 



(7.13) 



Expanding eqs.( |7.12] ) and ( |7.13| ) above as power series of r, we have access to the 9 expan- 
sion of Zj of eq.( [7.10|) by writing r = j9, namely 



+ 9' {}j' + ^f {t+l)+Oif)] + 
= 2 - ( I - ^(t + 2) + O (1) ) - (^If + Ijit + 2) + 



(7.14) 
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Prom the = (1 — gj Qci^l'^ terms, we read off the singular parts of Zj and dZj/dai 



In particular, the leading contribution to -^jUng does not depend on the a's and p's. To 
evaluate {nij) at large j, we now note that when all a's and p's are one, we have 

Z,({1;1})|.„., ^(m,) [ M^" (Q-) 

where we have used again the values of l})|sing above. The above result for quad- 

rangulations is in agreement with that of Ref.0 (eq.(4.19)). Using the above values of 
(mi), we therefore obtain by differentiating with respect to j: 



(7.17) 




at large j. Gathering eqs. (|7.8| ), ( [7.15| ) and ( [7.17| ), we finally deduce that dV /dai = dt/dai 
as announced. 

7.3. Applications 

As a first application of the above scheme, let us re-derive the results of Sections 4.1 
and 4.2 for /c = 1 with ai = a and pi = p. Eqs. (|3. 121 ) reduce to Zi{l — ga{Zi + Z2)) = pa 



(Q.) and Zi{l - gaZ2) = pa (E.T.). Substituting Zi = 2{t + l){t + 4)/((t + 2)(t + 3)), 
Z2 = 2(t + 2)(t + 5)/((t + 3)(t + 4)) and g = 1/12 (Q.) or Zi = 2(t + l)(t + 5)/((t + 2)(t + 4)), 
Z2 = 2{t + 2){t + 6)/{{t + 3){t + 5)) and ^ = 1/8 (E.T.), we get the equations for t: 

6{t + l){t + 2)(t + 4) - a{2{t + l){t^ + 6t + 6) + 3p{t + 2f{t + 3)) = (Q.) 

4(t + l){t + 3)(t + 5) - a{t + 2)((t + l){t + 6) + 2p{t + 3)(t + 4)) = (E.T.) 

(7.18) 

with to = —1 in both cases at a = 0. We finally substitute t = F — 1 in eqs. (|7.18|) to 



recover eqs. ( [4.9| ) and (|4.18| ). 
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To conclude this section, let us now display the algebraic equations as obtained via 
the above scheme for T in the case k = 2, with ai = a, a2 = 1, pi = p and p2 = cr- They 
read respectively for quadrangulations: 

6u^r{ur + 4){ur-u + 2){ur + u + 2){ur + 2u + 2) 

- a (u^T{uT + 4) (8 + Au{2 + 3T) + 2u^{3T^ + AT - 3) + (T^ + 2T^ + SF - 4) 

+ 2w^(r + l)(3r + 1) + u^T{T + 1)2) + 3p{uT + 2){uT-u + 2f{ur + 2u + 2)^ j = 

(7.19) 

with u = \JA — 3cr, and for Eulerian triangulations: 

i92r(r + 'u)(2r + u - i)(2r + w + 3) 

- a{2V + w - 3)(2r + w + 5) (^48r(r + w) + p{u^ - i)(2r + w + i)(2r + w - 1)) = 

(7.20) 

with u = a/(27^^^ct)7(3^^^^- Again these equations are supplemented by the condition 
that r = 1 when Q: = p = cr = lor alternatively that F = when o: = 0. 
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Fig. 15: The probabilities P{n) of having n nearest (resp. next nearest) 
neighboring vertices in infinite quadrangulation (a) and Eulerian triangula- 
tions (b), for n = 1,2,..., 20. We note a maximum probability reached in 
the case of next nearest neighbors at n = 3 for quadrangulations and n — 4 
for Eulerian triangulations, while the probabilities for nearest neighbors are 
strictly decreasing. 
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These equations allow to obtain various moments and probabilities involving next 
nearest neighboring vertices in infinite quadrangulations or Eulerian triangulations. For 
instance, we obtain the following first few moments of the numbers of next nearest neigh- 
boring vertices of the origin: 

, , 54 , 24948 , 17067888 

, , 36 , 2148 , 4668948 ^ , ^ ' 

The probabilities P{n) of having n next nearest neighboring vertices are displayed in 
Fig. 15 (a) for quadrangulations and (b) for Eulerian triangulations, and compared with 
the corresponding probabilities for nearest neighbors as given by eqs. (|5.2|) and ( |5.6| ). 



8. Conclusion 

A completely open question regards the generalization to other graphs than quad- 
rangulations and Eulerian triangulations. In Ref. 0, we have investigated a number of 
possible generalizations, either by allowing for faces with higher valences or by considering 
colored graphs, all of which turn out to have generating functions governed by discrete 
"integrable" recursion relations on the (suitably defined) geodesic distance, and for which 
compact "soliton-like" expressions also exist. These were derived by use of a somewhat 
different class of bijections between the duals of these graphs and some classes of so-called 
blossom trees with constrained valences [|TT|,[7[]. To imitate the approach of the present 
paper, one would first need to derive a bijection involving geodesic distances, and with ap- 
propriately labeled trees with arbitrary valences, a task not performed yet. Another more 
promising approach would consist in adapting the techniques of the present paper directly 
in the language of blossom trees. For quadrangulations and Eulerian triangulations, both 
described by binary blossom trees with buds and leaves (see Refs.|]Tl|,|T^ for details), this 
involves weighting leaves, vertices and edges with suitable weights pi (per leaf at the right 
of a face at depth i in the dual of the original graph), g (per vertex), and aj (per edge 
bordering a face of depth j in the original graph). Repeating this weighting procedure 
in more general cases of blossom trees should give access to the local environment in so- 
called constellations, for instance two-constellations, i.e. graphs with arbitrary even-sided 
faces, which generalize quadrangulations, or three-constellations, which generalize Eule- 
rian triangulations by allowing white faces with valences multiple of 3. Presumably, the 
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generalization of the blossom-tree construction proposed in Ref. [13| for arbitrary planar 



graphs should also allow to explore the local environment of a vertex as well. 
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